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Abstract 

We study grand unified theories based on an SU (5) x SU (5) gauge group in which the GUT 
scale, MguT) is the VEV of an exact or approximate modulus, and in which fast proton decay is 
avoided through a combination of a large triplet mass and small triplet couplings. These features 
are achieved by discrete symmetries. In many of our models, -Mgtjt i s generated naturally by the 
balance of higher dimension terms that lift the GUT modulus potential, and soft supersymmetry 
breaking masses. The theories often lead to interesting patterns of quark and lepton masses. We 
also discuss some distinctions between grand unified theories and string unification. 
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I. INTRODUCTION 



Lacking direct evidence, the unification of couplings is one of the few experimental hints 
both of supersymmetry and of unification of the gauge interactions of the standard model. 
The separation of the unification scale from the Planck scale has long suggested that the 
methods of effective field theory might be appropriate in understanding unification. Prob- 
ably the most troubling puzzle of such grand unified field theories (GUTs) is why Higgs 
doublets should be light, while their colored partners are massive. 

String theory offers an interesting perspective on these issues. In string theories, one 
often has unification of couplings, even though there is typically no scale at which one can 
speak of a four dimensional unifying gauge group. The unification scale, instead, is to be 
identified with the radius of some internal space, or with the string scale. In most cases, this 
scale is the threshold, not for a finite number of states, but for an infinite number of states. 
One of the appealing features of string models, stressed even in textbooks M, is that they 
can readily produce light doublets without light triplets. 

String theory also sheds light on some of the traditional questions of model building. In 
particular, it was long suspected that one should not expect continuous global symmetries 
in a quantum theory of gravity, and indeed there are no such symmetries in string theory 
@. String theory does frequently yield rather intricate patterns of discrete symmetries. 

This discussion suggests that we make a distinction between Grand Unification, i.e. the- 
ories where, for some range of energy scales, there are a finite number of fields and the 
gauge groups of the standard model are unified in a larger group, and String Unification, 
where the couplings are unified, but there is no such range of energies. There are a number 
of reasons to pursue Grand Unification, even if one imagines that the underlying theory is 
string theory. Perhaps the most important of these is coupling unification itself. Coupling 
unification is a robust property of the weakly coupled heterotic string. It is also true of the- 
ories in which some larger gauge symmetry is broken by Wilson lines (e.g. strongly coupled 
heterotic string, G 2 compactifications, etc.). But it does not hold more generally, and it is 
thus not clear in what sense it is a prediction of string theory. Certainly one could imagine 
that the explanation of coupling unification, even in string theory, is Grand Unification. 
Second, even though the appearance of massless Higgs doublets is an impressive feature of 
string theory, the absence of proton decay suggests that discrete symmetries may play an 
important role in the structure of the theory and, as Witten has recently stressed ||, such 
discrete symmetries might provide an alternative explanation of the presence of massless 
doublets. 

So even from the perspective of string theory, it is interesting to explore conventional 
grand unification, with particular emphasis on the role of discrete symmetries. From this 
perspective, however, there is a second puzzle: the origin of Mqut- In field theory models of 
unification, this scale often appears as an explicit input parameter. One might expect that 
in a string theory construction, the Higgs fields required to break the grand unified group 
would be massless in some approximation. In order that they obtain large expectation 
values, they would need flat or nearly flat potentials. This need not necessarily be the case 
- after all, a parameter of only 10 -2 is required in order to obtain this scale from the Planck 
scale. Still, the existence of a flat or nearly flat direction gives, as we will see, an appealing 
picture, allowing for a natural mechanism of generating the GUT scale in the effective field 
theory. This possibility will be the focus of much of this paper. 

We consider models in which Mqut is given by some combination of the Planck scale 
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and the supersymmetry breaking scale ||. The GUT group, Ggut, in our models, is broken 
by the expectation values of some fields, which correspond to an exact or approximate flat 
direction in the limit of unbroken supersymmetry. The location of the minimum, and the 
value of M G ut, are ultimately determined by supersymmetry breaking effects. 

It is easy to generate such an approximately flat direction by the use of discrete sym- 
metries. With symmetries such as a Z w , the potential along this direction arises from 
non-renormalizable terms suppressed by a high power of the Planck scale, such that, when 
balanced against a supersymmetry-breaking soft mass for the GUT breaking field, the cor- 
rect value of Mqut is obtained. 

Alternatively, Mgut may be an exact modulus of the theory. This can easily be achieved 
with a continuous global symmetry ||, but such global symmetries are unlikely to arise in 
a theory of gravity. It is possible to obtain exact flat directions with discrete symmetries, 
but typically one finds additional light fields which spoil the predictions of unification. We 
explain these difficulties, and construct some examples of exact flat directions which are 
almost, but not entirely, satisfactory. 

The discrete symmetries of our models also play a role in addressing the first puzzle 
mentioned above, namely, the presence of colored partners of the Higgs doublets. In a 
conventional grand unified model, because of the GUT symmetry, the Yukawa couplings of 
these partners are related by the gauge symmetry to standard-model Yukawa couplings of 
quarks and leptons, giving rise to fast proton decay. In the models we study, following ideas 
of H ||, [7|, a discrete symmetry distinguishes the triplets from the doublets. This symmetry 
not only allows a triplet mass while forbidding a doublet mass, but, as stressed in ||, may 
also suppress the triplet couplings to matter fields relative to the usual Yukawa couplings. 
Thus, proton decay is avoided by splitting not just the masses but also the couplings of the 
doublets and triplets. 

It is important to note that, as explained recently in ||, gauge coupling unification in 
minimal SU (5) implies an allowed window for the triplet masses, which is somewhat below 
Mgut- Proton decay, on the other hand, implies a lower bound, which lies above this window. 
The idea of Yukawa splitting between triplets and doublets may avoid this problem in more 
general models. (For another mechanism, see ref. 0.) 

As articulated by Witten, a symmetry that distinguishes doublets from triplets requires a 
semi-simple group, and indeed our models have SU(5) x SU(5) as Ggut- Still, the standard 
model gauge group lies in the diagonal SU(5), so that the main achievements of grand 
unification are maintained: hypercharge is quantized, and all standard model couplings 
originate from a single coupling, that of the diagonal SU (5). No large couplings are therefore 
required in order to ensure unification. In addition, the models require relatively few fields 
to achieve the initial stage of symmetry breaking, so that they remain weakly coupled up to 
the Planck scale. 

We review the mechanism of doublet-triplet splitting by discrete symmetries of |3|, || ffl] in 
Section II and describe the basic SU (5) x SU (5) theory and flat direction we consider. We 
explain how Mgut may arise naturally along an approximate flat direction in Section CI. 
Some explicit models that realize these two ideas are exhibited in Section [TV]. We argue that 
these are the simplest models which satisfy all of our requirements: along this direction, 
all these fields get Mgut masses, and a discrete symmetry, which forbids a doublet mass, 
remains unbroken. The models contain two pairs of bifundamentals, three adjoints and a 
gauge singlet. As we show in Appendix A, within such models a continuous symmetry is 
required to obtain an exact flat direction. 
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In Section [V|, we show that an exact flat direction, which is nearly satisfactory, can be 
obtained by considering a different pattern of expectation values. However, in these models, 
not all fields with standard model quantum numbers gain mass. In fact, a full SU(5) adjoint 
representation remains light, so that SU(3) becomes strongly coupled near M GU t- 

In the theories we consider, the standard-model Higgs doublets may originate from a 
5 and a 5 either of the same SU (5) factor, or of different SU (5) factors. To cancel SU (5) 
anomalies, one can then arrange the standard model matter fields in different representations 
of 577(5) x 577(5), or add an extra 5 + 5 pair. We consider these different possibilities in 
Section |V|. It should be emphasized that the discrete symmetry must be broken at some 
stage in order to allow a fi term for the standard model Higgs doublets, and we discuss some 
possibilities for accomplishing this. In particular, we show that a Higgs doublet mass may 
be generated automatically once the discrete symmetry is broken, without having to add 
any new Higgs couplings. 

Since, in this framework, the standard model matter and Higgs fields may originate from 
representations of the different S77(5)'s, some of their Yukawa couplings arise from non- 
renormalizable terms. This situation introduces a small parameter of order Mqut/^pi ~ 
10~ 2 into the Yukawa matrices. As we show in Section |VI]] (with further details given 
in Appendix 0), the resulting quark and lepton mass matrices are viable. Furthermore, 
some of the mass hierarchies are automatically explained by the gauge quantum numbers. 
When the standard Yukawa couplings of the Higgs doublets arise from non-renormalizable 
terms, the Yukawa couplings of of their colored partners can be highly suppressed. This 
Yukawa- splitting is also explained in Section |V11. 



The framework described above, of an effective theory in which Mqut corresponds to 
an exact or approximate modulus, as we have argued, seems a plausible outcome of string 
theory. String theory typically leaves many light, undetermined moduli, and possesses dis- 
crete symmetries. We discuss some distinctions between field theory unification and string 
unification in Section |V111| . 



II. SPLITTING TRIPLETS FROM DOUBLETS 

In this section and the next, we describe the key elements of our models. As explained in 
the introduction, we wish to distinguish between Higgs doublets and triplets by a discrete 
symmetry. (Here and in the following, we use "Higgs" to denote the standard model doublets 
and their triplet partners, as opposed to GUT breaking fields.) As discussed most clearly 
by Witten ||, if one wants to explain the masslessness of Higgs doublets through discrete 
symmetries, there should be an unbroken discrete symmetry (at least at a high energy scale) 
which acts differently on doublets and triplets, and is not a subgroup of hypercharge. This 
is not possible with a single 577(5), 5*0(10) or E$, but it is not difficult to achieve with the 
group 577 (5) x SU (5), and other direct product groups, provided one has fields transforming 
as bifundamentals. In this paper, we study the simplest such case, with an 577(5) x 577(5) 
gauge group. Here, one wants a symmetry which is a linear combination of a discrete 
symmetry acting on the bifundamentals, and a gauge transformation acting in one of the 
577(5) 's. The latter transformation can be taken to lie in the first 577(5), without loss of 



4 



generality (by combining with a hypercharge transformation), 
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Here a is an iV'th root of unity, and for simplicity we have taken N odd (for iV even it is a 
simple matter to modify our formulas). 



N+3 T 

a 2 $ 



2- 



Suppose that one has two pairs of bifundamentals, <3>j 
tential respects the Z N symmetry 

$i — > $i — > a~ 1( 5i, $ 2 

Then the vacuum expectation values (VEVs) 
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preserve a symmetry which is a combination of the original symmetry of eq. (0) 
and the hypercharge transformation of eq. (JI|). Because the symmetry is unbroken, this 
structure of expectation values is natural; it is automatically an extremum of the potential. 
We will discuss another possible set of expectation values, with only and ($2) nonzero, 
in section [V[ 

In fact, Barr constructed a model with precisely these features, for the case Z N = Z 2 [0. 
The superpotential in his model is: 

W = Mi($i<ii + $2*2) 

+ -!- [($!$!$!$!) + ($ 2 I>2$2$2) + ($1$1$2$2) + ($l$2$2$l) 

+ ($i$i)($i$i) + ($2$ 2 )($2$2) + ($l$l)($2$2) + ($l$ 2 )(*2*l) 

+ ($i$2$l5 , 2) + ($2$l$2$l) + ($1$2)($1$2) + (4 , 2$l)(4 , 2$l) • (4) 

(For simplicity, we have not distinguished the different couplings.) It is easy to see that 
for such a superpotential, there is a solution of the form (|3|). Moreover, one can also see 
that all fields gain mass. Note, in particular, that the only continuous symmetry of the 
superpotential is SU(5) x SU(5). As a result, if we turn on the gauge coupling, the only 
massless fields we expect to find are those associated with the Goldstone bosons for the 
breaking of this symmetry and their superpartners, which are Higgsed. 

This superpotential is not of the form we seek. The VEVs (§) are not a flat direction of 
the superpotential. Instead, the potential has a minimum with VEVs of order Mqut if M\ 
and M 2 are of order M GU t- ^gut appears as an explicit input parameter; both as mass 
terms (Mi), and as the scale suppressing higher- dimension terms (M 2 ). This superpotential 
is puzzling if the $'s are the only Mqut fields. Generating these couplings in a theory with 
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only fields of O(Mqut) -masses and marginal or relevant couplings requires the addition of, 
at least, a (24, 24) [[J with mass of order Mqut and couplings to the $ fields. 

However, for the purposes of splitting doublets and triplets, Barr's model is the simplest 
example of the mechanism we will use. It is now clear how triplets can gain mass while 
doublets remain massless. Assigning charges to the Higgses so that they cannot gain 
mass means that they can have couplings to $i or $ x , but not to $2 or $ 2 . From the point 
of view of the low energy theory, it is easy to choose Z N charges for the Higgses, such that 
the unbroken Z' N will forbid a doublet mass. 

If we want the unbroken Z' N to distinguish between doublet and triplet mass terms, the 
standard model Higgses must arise from fields charged under different SU(5) factors. 1 The 
standard model Higgses therefore come from, for example, 

h(5,l), h\l,5). (5) 

(Here and in the following, unprimed (primed) fields transform under the first (second) 
SU(5).) This possibility was discussed in 0. We refine this statement and consider addi- 



tional possibilities in Section VI 



The fact that the unbroken Z' N symmetry in the low energy theory below Mqtjt dis- 
tinguishes Higgs doublets and triplets may be used to relax the bound on the triplet mass 
from proton decay. Suppose that triplets are lighter than M GU t- Since triplets and doublets 
have different charges under the unbroken Z' N of the low-energy theory, their couplings to 
standard-model matter fields are typically not the same. From the point of view of the full 
theory, if, as mentioned above, the Higgs doublets come from a 5 and 5 of different S77(5)'s, 
then some standard-model Yukawa couplings originate from higher- dimension terms, involv- 
ing the bifundamentals. Thus, the triplets couplings to matter fields may be suppressed, 
so that their contribution to proton decay is small. Of course, the triplet mass has to be 
sufficiently high in order to preserve the unification of couplings. But in minimal 577(5) 
models, this requirement actually gives both an upper and a lower bound on the triplet 
mass. While in our models the triplet masses arise from renormalizable couplings to $1 
and $1, it may be possible to construct models in which these masses are lower, as low as 
permitted by gauge coupling unification. 



III. THE ORIGIN OF M GUT 

We are interested in models in which Mgut corresponds to an exact, or nearly exact flat 
direction. The main difficulty in constructing such models is that, while no potential should 
be generated along this direction (or at least no sizable potential), all fields charged under 
the standard model gauge group should acquire C(Mgut) masses. As we will see, however, 
it is possible to construct SU(5) x SU(5) models that are fairly simple yet possess these two 
features. 

We note that various constructions of this type have been put forward in the past. For 
example, a model based on 5*0(10), which also solves the doublet-triplet problem (using 
the so-called "Dimopoulos-Wilczek mechanism") was presented in || 2 . In this model, the 



1 For h(5, 1) and h(5, 1), the doublet and triplet mass terms /13/13 and h^hi have the same Z' N charges. 

2 Other examples have been given in jUj [yj [T^]. However, in most cases, these models were only studied 
at the level of renormalizable terms. 
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existence of flat directions, as well as the masslessness of the Higgs doublets, are insured by 
global continuous i?-symmetries, which restrict the Lagrangian to renormalizable terms. The 
pattern of expectation values which give massless doublets does not respect any particular 
symmetry; it is an accident of the restriction of the model to renormalizable terms. Allowing 
only discrete symmetries, the flat direction is no longer flat, and the desired pattern of 
expectation values requires fine tuning. It should be noted that, even with the use of global 
symmetries, the construction requires six adjoints, two symmetric tensors and two spinor 
representations just in order to obtain the first stage of symmetry breaking. 

If one does not worry about the spectrum, it is not difficult to obtain models in which 
the flat direction is an exact modulus using discrete R symmetries. For example, suppose 
one has a conventional £77(5) theory with a single adjoint. Requiring that the adjoint, A, 
be neutral under the R symmetry, while the superpotential transforms with a non-trivial 
phase, forbids any superpotential for A. So the potential for A is flat. However, this simple 
example also illustrates the main difficulty: if we choose to break to SU(3) x SU(2) x U(l), 
there are additional massless fields transforming as (8, 1) and (1,3) under SU(3) x SU(2). 
This difficulty arises also in the case of approximate flat directions. E.g., if one has an 
ordinary (non-i?) Ziq under which A — > e~^A, then the leading term in the superpotential 
involving A alone is A 10 , so the flat direction is very flat, but the octet and triplet are very 
light. 

We will consider, as above, product groups. Most of our models will involve two pairs of 
bifundamentals, as described in section [TV], with the VEVs (|3]). To avoid the problem of light 
fields, we will study models with additional fields beyond the minimal set of bifundamentals. 
The Zn global symmetry we discussed does not forbid a superpotential in the would-be flat 
directions. In fact, all nonzero gauge invariants, such as $i$i, $2^2, and $^$2 (contracted 
with e-tensors) are also singlets, and can appear in the superpotential to any power. 
It is easy to forbid such terms using an additional continuous symmetry in order to obtain 
an exact flat direction. Discrete i?-symmetries can give exact flat directions involving the 
$ fields but, as in our simple ££7(5) example, there are unwanted massless fields in these 
directions. As we will show, with the simplest possible field content, one cannot obtain 
exact flat directions with all fields massive using only discrete symmetries (continuous global 
symmetries can do the job). 

But the fact that the flat directions are only approximate can be a virtue, yielding a 
simple mechanism for generating Mqut 0- Suppose that the first contribution to the F 
terms comes from the superpotential term 

where X stands collectively for the fields of the theory and Mpi is the Planck scale. Suppose 
further that, once supersymmetry is broken, X acquires a negative soft mass squared, m 2 . 
Then the potential for X is 

V = ~m W + ^lXr 2 , (7) 



leading to an A-VEV of the order 



1 

Til \ n-2 



(X) ~ ( j^-J M P1 . (8) 
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For n around 10, this is around Mqut- 

We therefore aim, in our models, for superpotentials which generate F terms only from 
nonrenormalizable terms such as ($i$i) 5 , or (^f^) 2 . 

In models in which Mqut corresponds to an exactly flat direction, we will assume that 
Mqut is generated by supersymmetry breaking effects. Explicit examples of this type were 
studied in |T3 in the context of gauge mediated supersymmetry breaking, and in WA 



in 

the context of anomaly mediated supersymmtery breaking. The mechanism of [13 leads 
however to large flavor violation, unless extra structure is invoked. We note that other 
mechanisms for generating Mqut exist in the literature [II], [12]. In these models, Mqut is 
related to the dynamical scale of some gauge group. 

In the context of string theory, there is another mechanism which one might imagine 
could generate M GUT . This is the appearance of a Fayet-Iliopoulos term, with a coefficient 
of order , perhaps small enough to explain the scale ||15|| . In this case, Mqut might not 
be related to supersymmetry breaking. It is still necessary to have directions in field space 
which are flat or nearly flat with respect to the unified gauge group and the superpotential. 



IV. MODELS WITH 'MESONIC FLAT DIRECTIONS 

We now turn to some examples which demonstrate the various ideas we have discussed. 
As described in the introduction, our models contain two pairs of bifundamentals and 
i = 1,2. We impose a global symmetry, which ensures the masslessness of the 
Higgs doublets. The charges of the bifundamentals are listed in table |. Note that the 
bifundamentals form vectorlike representations of the Z^ - the charges of $i and $i sum 
to zero, and similarly the charges of $2 and $2- 

We would like to find a vanishing potential along the direction discussed in the introduc- 
tion [eq. (D]. Along this direction, SU(5) x SU(5) x Z N breaks into SU(3) D x SU(2) D x 
U(1)d x Z' N . Z' N is a combination of the original discrete symmetry and a discrete "hy- 
percharge" subgroup of SU (5)i. Various components of $j, $j are uneaten. In particular, 
the following standard-model SU(3) x SU{2) representations appear: 3 x (8, 1), 3 x (1,3), 
2 x (3,2) and 2 x (3,2). 

It is easy to see that we cannot give mass to the bifundamentals by couplings among the 
bifundamentals alone, without spoiling the flat directions. Adding two fields transforming 
as adjoints under the diagonal SU(5) forces Zn charges that are inconsistent with an exact 
low-energy Z' N and, furthermore, leaves light, incomplete SU(5) multiplets. At a minimum, 
we need to include in the model three fields transforming as adjoints under the diagonal 
SU(5). We add three adjoints of SU(5)i, A i=12 , 3 , as well as a gauge singlet S, with the 
superpotential 

W = Ai 2 $iAi<i> 2 + \ 21 <£> 2 A 2 $ 1 + An$iA 3 $i + A 2 2$2^35 , 2 + rjiaSAxAi + r] 33 SA 3 A 3 . (9) 

This superpotential indeed preserves a Z^ symmetry, for appropriate choices of the singlet 
and adjoint charges. The field content of the model is summarized in Table Q. The fields 
appearing below the first horizontal line contain the Higgs doublets. Most often, these will 
come from h and h', which we always include. The remaining two fields, h' and h, may or 



may not appear. We postpone further discussion of these Higgses until Section [VT], and for 
now concentrate on the GUT breaking fields. In Appendix B we show that the field content 
in Table | is still small enough that there are no Landau poles below the Planck scale. 



Field SU(S) x SU(5) x Z N 





(5,5,1) 




(5,5,-1) 




(5,5,(7V-3)/2) 




(5,5,(7V + 3)/2) 


Ai 


(24, 1, (iV - 5)/2) 




(24, 1, (iV + 5)/2) 


^3 


(24,1,0) 




(1,1,0) 




(5,1,1) 


b! 


(1,5,0) 


h 


(5,1,0) 


h! 


(1,5,-1) 



TABLE I: Basic fields and their charge assignments 

Indeed, the potential vanishes for 

($i> = ($!> = Vl x diag(l,l,l,0,0), 
($ 2 > = ($ 2 > = v 2 x diag(0,0,0,l,l), 

(5) = s, (10) 
(A) = 0, 

and all pseudo-Goldstone bosons get mass. The flatness condition requires 3 

Anfi = \22vl. (11) 

For superpotential couplings of order one, the two VEVs are not very different. 

Clearly, however, various terms which are allowed by the symmetries would spoil the 
flatness of the potential. For example, any power of $i$i, $i$2, or generates a potential. 
In the following, we consider different symmetries that forbid all or some of these terms. 



A. Exact flat directions 

If one allows continuous global symmetries, it is easy to obtain models in which the 
direction ([T0| ) is exactly flat. The superpotential @ preserves a U(1)r symmetry under 
which all adjoints A{ have charge 1, the singlet S has charge —1, and the bifundamentals $j, 
$j have charge zero. The superpotential has charge 1. For completeness we list the charges 
of all fields in Table O. 

The U(1)r symmetry forbids all the 'dangerous' non-renormalizable terms that would 
otherwise lift the flat direction. It does so by requiring that in each allowed term there 
is at least one adjoint. All terms with two adjoints or with a single Ai or A 2 are trivially 



3 To see that one needs to include Lagrnage multiplier terms ensuring the tracelessness of the adjoints. 
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Field SU{5) x SU(5) x Z N x U(1) R 



*1 


(5,5,1,0) 




(5,5,JV- 1,0) 


$ 2 


(5,5,(iV-3)/2,0) 


$ 2 


(5,5,(7V + 3)/2,0) 


Al 


(24,l,(iV-5)/2,l) 


A 2 


(24, 1, (iV + 5)/2, 1) 


A 3 


(24,1,0,1) 


5 


(1,1,0,-1) 



TABLE II: Exact flat direction 

harmless. The non-renormalizable terms with a single A3 modify the flatness condition (|i~l~l); 
the precise relation between v\ and v% is then changed, but not the fact that there is an 
exactly flat direction of the form (^). 

If, however, we consider just a discrete R symmetry, the flat direction is lifted. The 
problem is easy to understand. To obtain an exact flat direction for the $'s, we would like 
them to be neutral under some R symmetry. But this fixes the R transformation laws of the 
Aj's and, as a result, of S. In particular, S cannot be neutral. We show in Appendix A that 
there are always terms of the form S ,n ($$) m which transform properly under the symmetries. 
So in this class of models, it is easy to obtain an exact flat direction for the bifundamentals 
using only discrete R symmetries, but it is not possible to obtain an exact flat direction 
simultaneously for S. In the next section, we discuss approximate flat directions. 

B. Discrete Symmetry: Approximate flat directions 

As we just saw, we can get a model with an approximate flat direction by replacing the 
continuous R symmetry above by a discrete R symmetry, Z^. (The field content and charge 
assignments of this model are presented explicitly in Table |V|.) Then the terms 

—^S M -\ —L-sS 1 *- 1 ®^ z=l,2, (12) 

are allowed and lift the flat direction. 

In this case, the S flat direction is lifted by a superpotential term of dimension M — 1, 
while the $ flat direction (here $ stands for all bifundamentals) is lifted by superpotential 
terms of dimension M + 1. As discussed in Section |TIT], if the singlet and bifundamentals get 
negative masses-squared by supersymmetry-breaking, they will be stabilized near Mqut f° r 
M ~ 9,10. 

We thus obtain a model in which the GUT breaking VEVs are naturally obtained at the 
right scale, all fields associated with the GUT breaking obtain masses of order Mqut, and 
the low energy theory has an unbroken global Z' N symmetry, which can distinguish between 
Higgs doublets and triplets. 

It is straightforward to generate other models that share these features. For example, 
consider giving the $'s nonzero charges under the R symmetry. We take the R symmetry 
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Field 517(5) x SU(5) x Z N x 



*! 


(5,5,1,0) 


$1 


(5, 5, A -1,0) 


$2 


(5,5,(A-3)/2,3) 




(5,5,(A + 3)/2,8) 


^ 


(24,l,(A-5)/2,4) 




(24,l,(A + 5)/2,9) 




(24,1,0,1) 


5 


(1,1,0,10) 



TABLE III: Model B: Approximate flat direction 



to be Z^ with the charges given in Table [QI|. The symmetries allow, among others, the 
following terms: 

The flat direction is lifted by terms of dimensions 9 and 10 and consequently a reasonable 
value for M GUT is obtained. 



V. MODELS WITH 'BARYONIC FLAT DIRECTIONS 

As we saw in the previous section, it is impossible in the models considered so far to 
obtain an exact mesonic flat direction using discrete symmetries only. This can easily be 
done, however, if we consider the following direction: 

= v x diag(l,l,l,0,0), 
($ 2 > = vxdiag(0,0,0,l,l), 

= 0. (14) 

Again, these VEVs break the gauge symmetry to the standard-model gauge group, and 
preserve the Z' N global symmetry irrespectively of the $j charge assignments [the Zjv-charges 
of $i and $ 2 are still 1 and (N - 3) /2]. 

Apart from the possibility of obtaining exact flat directions, these directions are also 
promising for obtaining potentials which are flat to a very high degree. The only non- 
vanishing gauge invariant combinations are powers of the baryon operator If the first 
such term appearing in the superpotential is the baryon squared, it can naturally generate 
the scale Mqut once supersymmetry is broken. 

It is easy to see that discrete symmetries can guarantee an exactly vanishing potential in 
this case. For example, with no additional fields in the model, we can choose a Z^ symmetry 
under which the <3>'s have charge zero and the $'s have charge one. Then, assuming the 
superpotential carries i?-charge 2, the only terms allowed contain two $'s, and the potential 
vanishes. 

It is, however, very hard to generate masses for all GUT breaking fields in these models, 
or even for all fields but one (or more) sets in a complete SU(5) representation. In particular, 
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note that the uneaten fields in the $'s are an S77(3) octet and an 577(2) triplet, and these 
can only get mass from the 'baryon' operators $f and $i$2- 

In the following, we present three types of models. Each of these models demonstrates 
some nice features that can be achieved with 'baryonic' flat directions, but also suffers from 
some problems. Model A is closest in spirit to the models with 'mesonic' flat directions 
discussed in the previous section. It has an approximately flat potential but a light 24 of 
S77(5). With such a large representation, the S77(3) coupling blows up below M G ut, unless 
the fields in the 24 get masses of at least 1 TeV. Model B has no light fields beyond the 
MSSM fields, but the flat direction is lifted by a dimension-5 superpotential term. Model C 
has an exact flat direction but, again, has a light 24. Both models B and C have no adjoints, 
but in addition to the light 24 have incomplete multiplets at a scale of order 10 12 GeV which 
might spoil coupling unification. 



A. Model A: Approximate Flat Direction and a Light 24 

The field content of our first model appears in Table |TV|. We take the superpotential to 
carry i?-charge 1. The symmetry allows the superpotential 



Field SU{5) x SU(5) x Z N x Z { 



*1 


(5,5,1,1) 


$1 


(5, 5, A -1,0) 


$ 2 


(5,5,(A-3)/2,l) 


§ 2 


(5,5,(iV + 3)/2,0) 


.4 


(24,1,0,0) 


#1 


(1,24, (A + 5)/2,0) 


B 2 


(1,24, (iV-5)/2,0) 


Si 


(1,1,-5,0) 


S 2 


(1,1,5,0) 



TABLE IV: 'Baryon' model A 



W = 5 2 $i$2 + J Bi$2$i + + A$ 2 $2 + SiBf + S 2 Bl 

+ w^ l)2 + Wi (SlS2){ ^ 22)2 + " ' ' (15) 

Here we are interested in nonzero VEVs for Si, S 2 as well. We only show terms that either 
give mass to some fields, or generate a potential. We get an approximate flat direction, 
allowing for a natural generation of Mqut- Note, however, that for general N, the 
symmetry is spontaneously broken by the singlet VEVs, and one needs to carefully choose 
the Higgs charges, so that the doublets remain light. In addition, one needs to ensure that 
higher-dimension terms involving the bifundamentals and the singlets, that destabilize the 



VEV fll4|) , are sufficiently suppressed. Alternatively, if we insist on having an unbroken 
global symmetry, we can work with N = 2, with only $2, $2, and Bi odd under the Z 2 - The 
unpleasant feature of this model is that it leaves a light 24. 
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B. Model B: Approximate Flat Direction with No Light Fields 

We can write down a much simpler model with just the $'s and $'s. We impose a Z 2 x Z 2 
global symmetry. The first Z 2 is our original with N = 2, under which only $ 2 and $2 
are odd. Under the second Z 2 the two $'s are odd. The allowed superpotential is: 

W = — ($i$i$i$i + § 1 <&2&l®2 + $2^1 < ^2 <: &1 + $2 $2 "I* 2 $2 + §i&i§ 2 § 2 + §i§ 2 &2®l) 

+ + + + ■ (16) 

The dimension-4 terms give mass to all components of the $'s (these therefore get mass 
two orders of magnitude below Mgut). The first two dimension-5 terms give mass to the 
SU (2)-triplet and SU (3)-octet in the $'s (these therefore get mass four orders of magnitude 
below Mgut, a potential problem for coupling unification). Finally, the problematic term is 
the last one; it lifts the flat direction already at dimension-5, and is thus inappropriate for 
generating Mgut- In order that the model is acceptable, we need to identify M = Mgut in 
eq. flTfil) , that is, put Mqut by hand. 



C. Model C: Exactly Flat Direction 

Consider a model with the bifundamentals and singlets, where the $j fields carry a Z§ r 
charge +1 and all other fields ($, and Si) are Z^-neutral. The superpotential has Z^-charge 
+2. Thus all terms must include two $.j-fields and consequently the baryonic direction (|HD 
is exactly flat. Consider additional symmetries that allow only the following terms up to 
dimension-5: 

W = j^- (S , i$i$i$i$2 + S , 2 $2$l < I ) 2$2 + S , 3 $i$i$2$l) • (17) 

This superpotential leaves a light 24, while all other fields [which include (8, 1) + (1, 3) plus 
complete 5'C/(5)-multiplets] get their masses at a scale of order Mq UT /Mpj which may be 
problematic for coupling unification. 



VI. STANDARD MODEL HIGGSES 

So far, we have concentrated on the GUT breaking fields. In Section we found models 
in which Mgut corresponds to an exact or approximate flat direction. Along this direction, 
all fields associated with the GUT breaking are heavy. In addition, there is an unbroken 
Z' N global symmetry. We will now see in more detail how this symmetry splits the doublets 
and triplets. For concreteness, we focus on the first model described in Section |IV B| , but 
it is straightforward to repeat this discussion for different versions of the model - only the 
charges under the discrete i?-symmetry will change. For convenience, we list in Table |V| the 
relevant fields. 

Let us next write down the most general renormalizable superpotential that involves the 
h fields: 

W x = h& x h' + h'Qih. (18) 
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Field SU{5) x SU{5) x Z N x 



$! 


(5,5,1,0) 


*1 


(5,5,iV-l,0) 


<D>2 


(5,5,(iV-3)/2,0) 


$ 2 


(5,5,(A + 3)/2,0) 


Ai 


(24, 1, (iV - 5)/2, 1) 


A 2 


(24, 1, (N + 5)/2, 1) 


A 3 


(24,1,0,1) 


S 


(1,1,0,M-1) 


h 


(5,1,1,1) 


h! 


(1,5,0,0) 


h 


(5,1,0,0) 


ti 


(1,5,JV-1,1) 



TABLE V: The fields and symmetries of the model of Section IV Bj 



Since the h fields do not couple to the $2 and S fields, only the triplets acquire masses. 

It is useful to look at the charges of the doublets and triplets under the unbroken Z' N of 
the low energy theory. We list these in Table [VI|. Clearly, no doublet mass term is allowed 



Field 


z' N 


h 3 





h 2 


(iV + 5)/2 


h 


1 


h 2 


(iV-3)/2 


h's 


-1 


h' 2 


-1 


% 





h> 2 






TABLE VI: Doublet (h 2 ) and triplet (h 3 ) Z^-charges. 

in the low energy theory below Mqut- Therefore, if we start from four Higgses, the theory 
contains four light doublets. 

It is easy to see that this result always holds if we insist on an unbroken Z' N and masses 
of order Mgut for all triplets. In order to give mass to all triplets, we need to allow both 
terms in W\ of eq. flTB|). Since $! and $1 have opposite Z N charges, the Z N charges of the 
Higgses must satisfy 

Q(h) + Q(h') + Q(h') + Q(h) = . (19) 

As a result, the terms 

h$ 2 h' , h'<S> 2 h , (20) 

have the same Z^ charges: they are either both allowed, or both forbidden. We are therefore 
left with four light doublets. There are then three options a priori: 
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A: The theory does not contain h' and h. ££7(5) anomalies are cancelled by appro- 
priate choices of £77(5) x ££7(5) representations for the standard model matter fields. 
The standard model Higgses come from fields charged under different ££7(5) 's. Some 
standard model Yukawa couplings arise from non-renormalizable terms. Z' N can be 
broken by supersymmetry breaking effects to generate the \x term. 

B: The theory does contain h! and h, but these remain massless 4 0. Witten speculates 
that these could be the messengers of supersymmetry breaking. The heavy triplets 
are from h and hi . 

— Bl: The standard model Higgses come from fields charged under different 
££7(5) 's, say, h and h', so that, again, some standard model Yukawa couplings 
arise from non-renormalizable terms. 

— B2: The standard model Higgses come from fields charged under a single ££7(5), 
say, h and h. Then, standard model fields can all be charged under the same 
££7(5), and all Yukawa couplings are renormalizable. 



C: The theory does contain h! and h. All triplets gain mass through the couplings (|TI 
The Z' N is broken at a high scale, so that one doublet pair also gets mass around M G ut- 
It is possible to arrange for an acceptable fi term for the remaining two doublets, for 



example, through the mechanisms proposed in ||T6| or in 0. This is most easily done 
by adding a gauge- singlet, Sh, charged under the Z N , with a GUT scale VEV. Then, 
in order to have masses of order M GUT , the doublets that couple to Sh must be charged 
under the same ££7(5), and the higgs doublets are charged under the second ££7(5). 

In this case, one should also make sure that the pattern of VEVs we are considering, 
eqn. (^), is not destabilized, that is, that $i, $i do not get VEVs in the last two 
entries. It is easy to see that such VEVs can be very small. For example, if Sh has 
charge 1 under the Z^, then the operator £^ +5 ^ 2 £ M-1 $i$ 2 ( I?i ( I ) i is allowed, and can 
generate the danegrous VEV. However, it is very suppressed. In fact, such operators 
might even generate a fi term of precisely the right size, without having to add any 
new coupling for the Higgses. The doublet mass will simply arise from the first term 
in ([ID. 5 

In the next section, we will study the implications for quark and lepton mass matrices. 



VII. YUKAWA MATRICES AND YUKAWA SPLITTING 

We now turn to the Yukawa couplings for quarks and leptons. We will see that some 
of the possibilities mentioned above are excluded, and that in others, the ££7(5) x ££7(5) 
gauge symmetry automatically suppresses some Yukawa couplings. Moreover, in many cases, 
triplet-Higgses couplings to matter are naturally suppressed, demonstrating the mechanism 
of 'Yukawa splitting". 



4 It is easy to forbid the relevant mass terms by choosing appropriate charges for h! and h. 

5 With a Zi instead of a Zn, one does not even have to break the Z2 in order to generate a large mass for 
the extra triplets. This is the mechanism employed in . 
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The chiral part of our model contains the Higgs and matter fields. Taking into account the 
various possibilities for the Higgs fields discussed in the previous section and the requirement 
of anomaly cancellation, we have the following options: 

A: Standard model higgses from (5, 1)# + (1,5)h, no additional fundamentals, and the 
standard model fermion generations coming from one of the following sets of representations: 

1. 3x (5, l) + 2x (10,1) + (1,10) 

2. 2x [(5,1) + (1,10)] + (10,1) + (1,5) 

3. (5,1) + 2 x (1,5) + 3 x (1,10) 

Bl: Standard model higgses from (5, 1)# + (1, 5)#, additional fundamentals in (5, 1)# + 
(1,5)h, and the standard model fermion generations coming from one of the following sets 
of representations: 

1. 3x [(5,1) + (10,1)] 

2. 2x [(5,1) + (10,1)] + (1,5) + (1,10) 

3. (5, 1) + (10,1) + 2 x [(1,5) + (1,10)] 

4. 3x [(1,5) + (1,10)] 

B2, C: Standard model higgses from (5, l)#+(5, 1)#, additional fundamentals in (1, 5)#+ 
(1, 5)h, and the standard model fermion generations coming from one of the following sets 
of representations: 

1. 3x [(5,1) + (10,1)] 

2. 2x [(5,1) + (10,1)] + (1,5) + (1,10) 

3. (5, 1) + (10,1) + 2 x [(1,5) + (1,10)] 

4. 3x [(1,5) + (1,10)] 

Note that the only renormalizable Yukawa couplings involve 

(10, 1)(10, 1)(5, 1) H , (1, 10)(1, 5)(1, 5) H , (10, 1)(5, 1)(5, 1) H . (21) 

Non-renormalizable terms will involve one or two bifundamental fields. The resulting Yukawa 
couplings are suppressed therefore by one or two powers of e = Mqut/Mpi ~ 10~ 2 . This 
leads to several interesting consequences. 

First, a phenomenologically viable model must have a renormalizable top- Yukawa. This 
requirement excludes model (3) of class A and models (4) of classes B and C. 

Second, there is a single possibility to have all Yukawa couplings coming from renormal- 
izable terms, and that is model (1) of classes B2 and C. Here the gauge symmetry plays no 
role in the Yukawa hierarchy. 

Third, when down-type Yukawas invlove matter fields charged under different 577(5) 's, 
say, a (10, 1) and a (1,5), they give either a down-quark or a charged-lepton mass term. 
This allows for the exciting possibility that, just by the gauge quantum numbers, bottom- 
tau unification is maintained, but similar relations for the two light generations are not! 
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Furthermore, we see that it is possible to explain some features of the Yukawa hierarchy 
by the SU(5) x SU(5i) symmetry. Take, for example, model (1) of class Bl. The mass 
matrices here are of the form 




M u ~ ((5,1) H ) 1 1 1 , M d ~ ((1,5V) 



fe e 

e e e | . (22) 

Thus, the small ratio m^/rrit is explained by the gauge symmetry. The hierarchy within 
each sector requires, however, some additional flavor physics. We show the mass matrices 
obtained in the remaining cases in Appendix |y. 

Finally, we can see how doublet-triplet Yukawa- splitting works. Take again model (1) of 
class Bl. The entries in of eq. (|22"D come from the non-renormalizable terms 



(10,l)(5,l)(l,5) H (5,5) i2 . (23) 

Note that these terms induce Yukawa couplings for the doublet Higgs fields only. Yukawa 
couplings for the triplet Higgs fields could have been induced by the following terms: 

(10, 1)(5, 1)(1, 5)^(5, 5) #1 . (24) 



However, if the discrete charge assignments are such that the terms in fl23|) are allowed, 
those in (E3) are, in general, forbidden. This situation can be understood also in terms of 
the unbroken low-energy Z' N symmetry that can allow Yukawa couplings for the doublets 
and forbid them for the triplets. Thus, the triplet- Yukawa couplings will only appear with 
Z^-breaking and can therefore be very strongly suppressed. As explained in section II, 
this mechanism of doublet-triplet Yukawa splitting is very helpful in solving the problem of 
proton decay. 



VIII. CONCLUSIONS: DISTINCTIONS BETWEEN GRAND UNIFICATION 
AND STRING UNIFICATION 

Various puzzles of grand unification find attractive solutions in models that combine 
product gauge groups and discrete symmetries. In this work we focused on a framework 
that combines an SU (5) x SU (5) gauge group with Abelian (Zfj) discrete symmetries. 

The gauge group breaks directly, at a scale Mqut, into the standard model SU(3) x 
SU(2) x U(l) group which resides in the diagonal £77(5). Let us first emphasize again 
that, in spite of the fact that there is no range of energies where there is a simple unifying 
gauge group, the attractive features of GUTs are maintained: there is coupling unification 
and hypercharge is quantized. At the same time, some of the problems associated with 
GUTs find solutions for which the combination of product groups and discrete symmetries 
is essential: 

(i) The doublet-triplet splitting problem has a symmetry-based solution. A discrete 
symmetry (that remains unbroken below Mgut) can forbid doublet masses while allowing 
triplet masses. At the same time, this symmetry can allow Yukawa couplings for the doublet 
while forbidding them for the triplet. Thus the inconsistency between the constraints from 
coupling unification and from proton decay that arises in minimal S77(5) can be avoided in 
this framework. 
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(ii) The GUT scale corresponds to an exact or an approximate flat direction. In the latter 
case, the value of Mgut is not an input parameter but is induced by an interplay between 
the Planck scale and the supersymmetry breaking scale. In addition to its intrinsic appeal, 
we argued that if the underlying theory is string or M theory, it is almost inevitable that 
M G ut should be a modulus. 

The combination of product gauge groups and discrete symmetries might have interesting 
implications for a third puzzle: 

(iii) A non-trivial flavor structure arises since, in general, some of the Yukawa couplings 
arise from non-renormalizable terms. Such couplings are suppressed therefore by powers of 
M G ut/Mpi. 

Our models demonstrate that it is possible to have fully consistent models of grand 
unification, with no light, extra incomplete (or even complete) SU(5) multiplets and with a 
relatively simple particle content at Mqut ; and in which the value of Mgut is explained in 
terms of the Planck scale and the supersymmetry breaking scale. 

The assumption that Mgut is a modulus opens up the possibility of a cosmological moduli 
problem. In the context of string theory, the difficulty is that moduli typically dominate the 
energy density at nucleosynthesis. Unless they are very massive, their subsequent decays 
reheat the universe only to temperatures of order a few keV, and the successful predictions of 



big bang nucleosynthesis are spoiled [|T7|, [Tj|. In the present case, however, the couplings of 
the modulus may be barely strong enough to evade this problem. Explicitly, the coefficient 
c in the leading contribution to the decay amplitude, Cm = ^m^ t -MFF ', can be 0(50) 
due to the large number of fields. 

In the introduction, we stressed that while string theory does not admit continuous 
global symmetries, it often yields discrete symmetries. These symmetries are usually gauge 
symmetries, and are subject to various consistency conditions. In weakly coupled string 
models as well as a range of strongly coupled string theories 0, [L9|, [20], |21| , it is possible to 
cancel all anomalies in discrete symmetries by postulating a discrete transformation law for 
some moduli fields. A priori, however, it is not clear that this provides a constraint on field 
theory models, since one might postulate the existence of several such moduli with Planck 
scale couplings. It is true that in most instances which have been studied, a single field 



can cancel all anomalies, but it seems unlikely that this holds in general f21 |. We have not 



explored possible anomaly constraints in our work, but this is certainly an issue worthy of 
further study. 

Whether or not of stringy origin, unification in the framework of product groups opens 
up new possibilities regarding the flavor structure of quarks and leptons and the question 
of proton decay. We have presented some aspects of the quark Yukawa hierarchy and of 
the doublet-triplet Yukawa splitting in Section [VII] . We are currently pursuing further 
investigations of these intriguing possibilities. 
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APPENDIX A: EXACT FLAT DIRECTIONS 



In this appendix we show that it is impossible to forbid all higher-dimension operators 
using a x discrete symmetry in a model with bifundamentals, adjoints and singlets. 

We denote the Z^-charge of each field by the corresponding small letter. For example, 
the charge of $j is <pi. The superpotential has charge to. All the equations involving the 
Zjd-charges are mod M, and all charges are integer. 

In order to allow the terms SA\ and $iA 3 <li of eq. (Q), the charges should satisfy 

4>i + 0i + a 3 = w , s + 2a 3 = w . (Al) 

We will now show that we can always find n, m so that 

1 S rt ($i$i) m (A2) 

is allowed. For this term to appear we need: 

n(w — 2a 3 ) + m(w — a 3 ) = w ==>■ (n + m — l)w = (2n + m)a 3 . (A3) 

This condition is satisfied for 

n + m-1 = qM , 2n + m = q'M , (A4) 

where q, q' are some integers. Now consider the charges under the Z^ symmetry. We can 
repeat the above, with w — > and M — > N. Therefore we get one additional equation for 
n, m: 

2n + m = q"N . (A5) 
The three equations (|A4|) -([A"5|) can be solved with q' = q" = kN, with integer k, giving: 

n = kNM - qM - 1 , m = (2q - kN)M + 2 . (A6) 

Since the charges are only defined mod(MA) we can choose: 

n = MN - qM - 1 , m = 2qM + 2 . (A7) 

A simple possibility is q = 0: the term 

S^- 1 ^,) 2 (A8) 

is always allowed. 

Note that if we have effectively just the Z^ symmetry, (as in the models of section [IV], 
where both S and are neutral under the Zn), then 

S* 1 - 1 ($!$!) 2 (A9) 

is allowed (corresponding to N = 1 in 
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APPENDIX B: LANDAU POLES 



Since our models require a large number of representations at Mqut, there is a potential 
problem of hitting a Landau pole below the Planck scale. We investigate this issue here. 
The RGE reads 

1 1 b M 2 

a{Mi) ~ a(M^j + 2^ ^Wi ( l > 

We will hit a Landau pole, a — > oo, at a scale M LP given by 



Mi 



LP 



M -^fe)' < B2 » 

Using «(Mgut) ~ 0.04 and requiring that the Landau pole is not reached below the Planck 
scale, that is M LP /M GUT > 10 2 , we find b < 34. 
For a supersymmetric SU(N) group, we have 

b = £C(0 a )-3C(adj), 

a 

C(fun) = 1/2, C(adj) = N, C(as) = N/2 - 1, (B3) 

where fun, adj and as stand for, respectively, the fundamental, adjoint and antisymmetric 
representations. Defining for an SU(5) x SU(5)' theory 

n h = #(5,1) + #(5,1), 
n 4 = #(5, 5) + #(5, 5), 
7U = #(24,1), 

n T = #(10,1) + #(10,1), (B4) 

we obtain for b(SU(5)): 

15 3 
b = -n h + + 5tia + -n T - 15. (B5) 

In particular, the maximal 6 arises if we all the MSSM fields (n^ = 5, ra^ = 3) and all 
M gut fields (n$ = 4, = 3) transform under the same SU(5): 

rih = 5, n$ = 4, = 3, rix = 3 ^> 6 = 17. (B6) 

Putting that in eq. ( |B2| ), with «gut — 1/25, we find 

M LP /50tt\ 4 

~ exp ~ 10 . (B7) 



M GUT V 17 

Thus our models are safe against Landau poles. 



APPENDIX C: QUARK AND LEPTON MASSES 

In this Appendix, we give expressions for the Yukawa matrices that arise in the various 
models we consider. The representations for the quark (or lepton) fields are given explicitly. 
The Higgs fields are denoted by h(5, 1), h{5, 1) and h'(l, 5). Bifundamental fields are denoted 
by $(5, 5) and $(5, 5). 
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For the up sector, the following combinations are relevant for inducing the Yukawa cou- 
plings: 



(10,1) (10,1) h, 
(1,10) (l,10)/i$, 
(10,1) (1,10) h $ 



(CI) 



For the down and charged lepton sectors, the following combinations are relevant for inducing 
the Yukawa couplings from h': 



(1,10) (1,5) h', 

(10,1) (5,1) h'Q, 

(1,10) (5,1)//$ [d/£], 

(10,1) (1,5) h! $ $ [d/4 (C2) 



and from h: 



(10,1) (5,1) h, 

(1,10) (1,5) 

(10,1) (1,5) h$ [d/£], 

(1,10) (5,1) h $ $ [d/4 



(C3) 



Entries marked with [d, £\ mean that these terms can give masses to either the down sector 
or the charged lepton sector but not to both. In the full high energy theory this depends 
on whether the bifundamental representations are $i or $ 2 (or $ x or $ 2 ). From the point 
of view of the low energy theory, below M G ut, this depends on the different Z' N charges 
carried by d^dn and £r£l- (^-breaking will eventually lift the zero masses.) 

Defining e = Mqut/Mpi, we get the following paramatric suppression due to the SU(5) x 
SU(5) gauge symmetry (entries in parenthesis vanish in either Ma or Mi): 

Model A(l), with (1, 10) as first generation: 




(C4) 



Model A(2), with (10, 1) as third generation and (5, 1) as first generation 



u 




(C5) 



Model Bl(l): 



(1 1 

M u ~ (h) 1 1 

\1 1 




(C6) 
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Model Bl(2), with (1, 10) + (1,5) as first generation: 



M u ~ (h) 



e e 2 e 2 ' 



1 1 
1 1 



M d ~ (^') 



/ 1 



Model Bl(3), with (10, 1) + (5, 1) as third generation: 

.2\ / 1 



M u ~ (h) 



e 2 e 2 1 



M d ~ (hf) 



Model B2/C(l): 



M„ 




M d ~ (/*) 



Model B2/C(2), with (1, 10) + (1,5) as first generation: 



/e e 2 e 2 



e z 1 1 
e 2 1 1 



M, 



'id 



(h) 



/e e e 2 
e e e 2 
\e 2 e 2 1 



/ e 





(e)\ 


e 




e 


J 


1 


(e)\ 


1 


(e) 


e 2 ) 


e / 




Model B2/C(3), with (10, 1) + (5, 1) as third generation: 





( f 2 


1 


1 


1 


1 


e 


(e 2 



(e) 




(C7) 



(C8) 



(C9) 



(CIO) 



(Gil) 



Note that, in principle, the Zjy symmetry may constrain these matrices further. We 
checked that none of the above models is excluded by these constraints. 
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